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Abstract 

In this paper we investigate the role of inelastic collision in the kinetic 
evolution of a highly overpopulated system starting from Glasma-type 
initial condition. Using the Gunion-Bertsch formula we derive the inelastic 
collision kernel under the collinear and small angle approximations. With 
both numerics and analytic analysis, we show that the inelastic process 
has two effects: globally changing (mostly reducing) the total particle 
number, while locally at small p always filling up the infrared regime 
extremely quickly. This latter effect is found to significantly speed up the 
emergence of local thermal form with vanishing local "chemical potential" 
and catalyzes the onset of gluon condensation to occur faster (as compared 
with the purely elastic case) in the overpopulated Glasma. 



1 Introduction 



Thermalization of the quark-gluon plasma is one of the most challenging prob- 
lems in current heavy ion physics. Starting with two colliding nuclei in a form 
of color glass condensate with high gluon occupation / ~ l/a s below satura- 
tion scale Q a [1, 2, 3] and following the initial impact, a subsequent strong field 
evolution stage (likely with instabilities [4]) till about the time 1/Q S is then suc- 
ceeded by a far-from-equilibrium gluon-dominant matter, the Glasma [5]. The 
evolution this Glasma stage toward a locally equilibrated quark-gluon plasma 
(QGP) is strongly indicated by phenomenology to be reached on the order of a 
fermi over c time (see e.g. [6]). Precisely how this occurs remains to be fully 
understood. Describing the pre-equilibrium evolution with kinetic equations is a 
very useful approach, based on which the so-called "bottom-up" therma lization 
scenario was developed [7, 8, 9] . There is however the complication of instability 
driven by anisotropy that may change this picture (see e.g. [10, 11, 12, 13]). 
There are also other kinetic-based approaches, see e.g. [14]. 

More recently an alternative thermalization scenario, based on crucial role 
of high initial overpopulation in the Glasma and kinetic evolution dominated 
by elastic collisions, has been proposed in [15, 16]. In this scenario, while 
the initial scale Q s is large compared with Aqcd and thus the coupling a s 
is small, the high occupation / ~ l/a s elevates the elastic scattering rate to 
be of the order O(l) rather than the usual 0{o? s ), and the Glasma is essen- 
tially an emergent strongly interacting matter with weak coupling albeit large 
aggregate of constituents. Two important scales are introduced to characterize 
the distribution, the hard cut-off scale A beyond which / << and the soft 
high-occupation scale A s below which / ~ l/a s . While the initial Glasma has 
the two scales overlapping A ~ A s ~ Q s , during the course of thermalization 
the two scales are separated eventually toward A s ~ a s A upon thermaliza- 
tion. One particularly nontrivial observation in the elastic-dominant picture 
is that the high initial overpopulation n/e 3 / 4 ~ l/a 1 / 4 » 1 and the con- 
servation of both energy and particle number will necessarily require the for- 
mation of a Bose-Einstein condensate that absorbs the excess gluons. This 
has been explicitly shown to occur by numerically solving the elastic kinetic 
equation derived under small angle approximation [16]. There have been in- 
tensive discussions related to this picture from a variety of approaches, see e.g. 
[17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. Strong evidences 
for the formation of such a Bose condensate have been reported for similar ther- 
malization problem in the classical-statistical lattice simulation of scalar field 
theory [19, 20, 21]. The case for non-Abelian gauge theory is more complicated 
and still under investigation [22, 23, 24, 25, 26]. 

One important question that has not been addressed in the above scenario is 
the role of inelastic processes. This issue could indeed be critical for at least two 
reasons (see discussions in e.g. [15, 16, 17, 18, 31, 32]). First of all the inelastic 
processes will spoil the particle number conservation, and one might naively 
argue that the excessive gluons in the ovcrpopulated Glasma could simply be 
eliminated by very fast inelastic collisions. Secondly, to make it even worse, the 
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inelastic processes are parametrically at the same order as the elastic processes 
(as opposed to naive power counting), so there appears no apparent dominance 
of the elastic over the inelastic and one may indeed worry that the inelastic 
could efficiently reduce total particle number. In such a situation, an explicit 
evaluation including both elastic and inelastic collisions becomes mandatory to 
clarify what will happen after including both t ypes of collisions. To be precise, 
once the inelastic processes are included, one does not expect any condensation 
in the ultimate thermal equilibrium because with long enough time the inelastic 
processes will always remove any excessive particles. The real question is, while 
in the purely elastic case the overpopulatcd system is driven toward a dynamical 
onset of condensation in a hnite time, whether the inelastic processes will even 
remove this transient off-equilibrium condensation or what else could happen. 

In this paper, we aim to address this important question by studying the 
kinetic evolution of a highly overpopulatcd system starting from Glasma-type 
initial condition with both 2h2 and 2f>3 scatterings. In Section 2 we will 
derive the inelastic collision kernel under the collinear and small angle approxi- 
mations using the Gunion-Bcrtsch formula. In Section 3 we will use numerical 
solutions as well as analytic analysis to understand the role of the inelastic pro- 
cess for both the global particle number change and the local behavior at small 
momentum region. Finally we will conclude in Section 4. As a first step toward 
understanding the inelastic contributions and for simplicity and unambiguity, 
we will focus on the static box case with isotropic distribution in this work and 
leave the study of expanding case for future work. 

2 Kinetic Evolution with Elastic and Inelastic 
Collisions 

In this section we will derive the kinetic evolution equation with both elastic 
and inelastic collisions. The kinetic equation deals with the gluon distribution 
function defined as: 

ft* \ (2?r) 3 dN 

The N g — 2(N% — 1) denotes the spin and color degeneracy factor. The Boltz- 
mann equation for fit, x, p) reads 

■D t f p =C 2 „2[f P }+C 2 „3[f P ], (2) 
where we denote /(i,x, p) by f p and 

V t = ^-d fi = d t + v p - V x (3) 
hp 

with v p = p/Ep and E p = |p|. For later convenience, we also introduce the 
following notations: 

9p = 1 + f P i h P = f P 9p = /p(1 + fp) ■ ( 4 ) 
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In what follows we will separately discuss the elastic term C202 and the inelastic 
term C 2 «3- 

2.1 The 2^2 process 

The collision kernel from the 2 o 2 process with full nonlinearity has been 
studied in [15, 16]. Here we only briefly summarize the main results. We have 
the 2—^2 collision matrix element square to be: 

|AW = 8/A^( 3 -5-f -£). (5) 

The 2f>2 collision kernel is then 

<W./p] = ~f T ±r\M 12 - >3p \ 2 (27r) 4 6 4 (p 1 +p 2 -p 3 -p) 

z J 123 ltj v 

x[(l + /„)(! + h)hh - / P /a(l + + /a)], (6) 



where 



The pre- factor 1/2 is a symmetry factor counterweighing the permutation of 1 
and 2 while the pre- factor 1/N g cancels the summation over the spin and color 
of gluon "p" in the matrix element square (5). 

The dominant contribution of 2 -H> 2 scattering in (5) comes from very small 
exchange momentum in t — > or u — > kinematic regimes, for which the 
incoming momenta only gets "deflected" by very small angle. If one uses this 
small angle approximation, then a rather neat kernel can be derived [16] : 



C 2 « 2 = ^A 2 s A^d p {p 2 



df p a s 
dp + A s Ip 



(8) 



1 

2 

with £ = J ^ coming from the leading log contribution. The hard scale A 
and soft scale A s in the above are defined via global integrals: 

/•OO 

AA 2 s /a 2 s = / dpp 2 f p (l + f p ) = I a , (9) 
Jo 

/•OO 

AA s /a s = / dpp 2 (2f p /p) = I b . (10) 
Jo 

For later convenience we also introduce the Debye scale defined as [33, 34] 

/■OO 

m 2 D = -a. / d P p 2 df p /dp = AA S . (11) 
Jo 

It is interesting to notice that in a weakly coupled thermal QGP one has the 
well-defined separation of scales, A ~ T, mo ~ gT, A s ~ <? 2 T. The matter 
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becomes strongly interacting when the scales "collapse" together. One way for 
that to happen is to have the system become really strongly coupled g — > 1 
which likely will be accompanied by change of underlying degrees of freedom 
[35] . The other possibility, as in the case of Glasma, is when the system is highly 
off-equilibrium and overpopulatcd / ~ 1/g 2 — in this case all the scales also 
become of the same order A ~ run ~ A s ~ Q s and make the system emerge as 
a strongly interacting matter. 

Clearly, both the full form in (6) and the small angle approximation 

form in (8) conserve the energy as well as particle number, as they should. In 
addition the Bose-Einstein distribution Jbe — l/[e( p_ ' i " T — 1] with any T and 
H (in correspondence to the two conserved quantities) is the fixed point solution 
that makes both (6) and (8) vanish. 

2.2 The 2^3 process 

We now turn to the collision kernel from the 2f>3 process as depicted in Fig.l. 
We denote the particle we are watching with momentum p, the softest external 
momentum with k, the exchange internal four- momentum with q^, and then 
the rest external momenta with Pi,2,3- The 2 •<-> 3 collision kernel can then 
be split into two pieces in which the particle p is on the two-particle side or 
three-particle side respectively (see Fig.l): 

C^ 3 = ~ / ^r\M lp « 23k \ 2 (2ir) 4 5 4 (p+ Pl -p 2 -p 3 -k) 

iV ff D Jl23fc 2 ^P 

x[(l + f p )(l + /O/a/a/fc - / p /i(l + h)(l + /s)(l + /*)], 
C b 2^3 = f ^r\ M ^ik P \ 2 (^) i S i (p+p 1 + k~p2^ P3 ) 

iV 5 4 Jl23fe ltj P 

x[(l + / p )(l + + f k )f 3 f 2 - f p hf k (l + / 3 )(1 + / 2 )],(13) 

where the gluon labeled by k will be treated as the soft emitted or absorbed 
gluon. The factor 1/6 counteracts the 6 equivalent permutations in 23fc in 
process l+p<H-2 + 3 + fc (see Fig. 1 left panel) and the factor 1/4 counteracts 
the 4 equivalent permutations in lk and 23 in process 2 + 3^l + k+ p (see 
Fig. 1 right panel). We note that the graphs in Fig. 1 are used to make the 
kinematics clear and it does not mean that only these two diagrams contribute: 
there are actually 25 different diagrams for C203 anc ^ 25 diagrams for C\^%. 
So the full matrix element square \Mi p ^2 3 k\ 2 is obtained by calculating 25 
Fcynman diagrams and it contains 6 equivalent kinematic setups in accordance 
with 6 permutations in 23fc (see Appendix A). We can then choose the kinematic 
setup corresponding to the Fig. 1 left panel, and multiply 6 to account other 5 
kinematic setups. Similarly, we can fix the kinematics for \M2 3 ^ikp\ 2 as in the 
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Figure 1: (Left) A typical Feynman diagram contributing to Cf-n^ P + 1 ^ 2 + 
3+fc. (Right) A typical Feynman diagram contributing to C\^ 3 : 1+k+p <-> 2+3. 

Fig. 1 right panel, and multiply 4 to get C\^ z . Thus we get: 



pa 



■'2<->3 



P3 - fc) 



^ / ^r|M 1 ^ 23fe | 2 (2 7 r) 4 J 4 (p+ Pl - p 2 

iV S Jl23fe 

X[(l + f p )(l + ./; ■h.hjl - f p fx(l + / 2 )(1 + / 3 )(1 + f k )\, 

W I ^r\M b 23 ^ lkp \ 2 (2n)^(p+ Pl + k-p 2 -p 3 ) 

ly g J 123k A£j P 

x [(1 + /„)(! + + f k )f 2 h - Uhfk{l + f 2 )(l + h)\, (14) 



where |M a,b | 2 are the matrix element square with the kinematics fixed according 
to Fig. 1. While the exact 2 o 3 matrix element square is known, it is very 
hard to be directly used in a kinetic approach. Following many previous studies 
involving this process [36, 37, 38, 39], we will use the so-called Gunion-Bertsch 
formula which is the collinear approximation and small angle approximation 
form of the exact matrix element and has been shown to give the dominant 
contribution in many cases. Leaving the technical details to the Appendix A, 
we quote the Gunion-Bertsch matrix clement square which is at the leading 
order in soft q and k expansion: 



l-^lp-(->23fe| 2 



|M 



23<->lfcpl 



64g 6 N^N : 
64g 6 N*N, 



(p-Pi) 



9 q 2 {q-k) 2 {p-k){ Pl -ky 
9 q 2 (q + k)2(p 2 -k)(p 3 -k)- 



(15) 



Like in the 2f)2 case, the matrix elements are dominated by the regime of 
very soft q and k. We can thus further simplify the collision kernel using similar 
small angle approximation as in the elastic case. There is though additional 
subtlety as now there are two soft scales. In fact, as shown by the detailed 
analysis in Appendix B, the whole collision kernel can be separated into two 
pieces corresponding to contributions from different kinetic domains. (Note 
that both graphs in Fig. 1 contribute to each of these domains.) In the domain 
with k being the softest scale i.e. the ultra-soft emission and absorption, the 
2 -<->• 3 essentially reduces to an effective 2 o 2 scattering with a collinear 
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splitting/merging, and the resulting contribution to the collision kernel becomes: 
-L f -^ r \M lp ^(2n) i 6 i (p + Pi-P2-l) 

iV 9 Jl2l ltj p 

x (g P 9ihfi - f P fi929i) V{\q = pi - p 2 |), (16) 



where the momentum labels I and q are as shown in Fig. 1, and we also intro- 
duced the term T>(\q = pi — P2I) arising from the splitting function integration. 
For the above form, one can further simplify the effective 2 <R- 2 part using small 
angle approximation as done in the clastic case. All the details are presented in 
Appendix B. On the other hand in the domain with q being the softest scale, the 
2 o 3 process effectively becomes a small angle 2f>2 scattering that brings one 
incoming particle slightly off-shell which undergoes a 1 -f-> 2 emission/combining. 
This part contributes the following to the collision kernel: 



—j Q dPlPlhl j_— xj 

dz 



dq 



~\ [9pf(i-z)pfzp — fp9(i-z)p9zp] 

[9p9zp/{l-z)fp/{l-z) - fpfzp/(l-z)9p/(l-z))] K17) 



(1-zY 



Note in the z-intcgration we introduce the upper cut z < z c : physically this 
is because the C^ 3 is from the domain with k being the softest external mo- 
mentum so the k < p condition would require z c < 1/2. In both kernels above, 
a number of infrared divergences appear. We will use the Debye scale mu as 
the infrared cut, e.g. J dq/q 3 ps ^/m 2 D = 1/(AA S ). We will also treat all lead- 
ing logs as order one constant. Again more detailed discussions are included 
in the Appendix B. Lastly, one can show that the inelastic kernel conserves 
energy while not particle number, and the fixed point solution (i.e. the equilib- 
rium distribution) is the Bose-Einstein distribution without chemical potential 
fBE = l/(e p/T - 1). 



2.3 The final kinetic equation 

Finally we combine the C202 and C2<->3 kernels, and the final kinetic equation 
under small angle approximation and collinear approximation reads 



where 



The expression for C%z+ 2 1S 



r cff 



T^tfp — £-2<->2[/p] + ^l<->2[/p]i 

°l<->2 — L '2«->3- 

dfp | h 
Bp I a 



Celt n 
2<->2 — <-"2-o-2 



L -2<-i-3' 



(18) 
(19) 

(20) 
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with D ~ 0(1) parameterizing the contribution from C^ 3 to effective 2 <-> 2. 
The expression for CJ^ 2 i s 

^ifU = {, a lR-r\ [ — [9pf(i-z)pfzp - f P 9(i-z)p9zp] 

+ J Q _ Z )4 Z [5pfl Z p/(l-*)/p/(l-*) - /p/zp/(l-*)Sp/(l-*))] j(21) 

where the constant i? ~ O(l) parameterizes the relative ratio of the order one 
constants between the elastic and inelastic kernels and the cut z c in z-integration 
should be small to be consistent with the kinematics k < p. The expressions for 
D and R can be found in Appendix B. 
A few comments are in order here: 

1) The kernel (20) conserves energy and particle number with fixed point solu- 
tion /be = l/[e( p_AI )/ T — 1], while the kernel (21) only conserves energy with 
fixed point solution Jbe — l/[e p / T — 1], so the total kernel (18) conserves only 
energy and the equilibrium solution should be Jbe = l/[e p/>T — 1] without any 
chemical potential which is different from the pure elastic case; 

2) In the near equilibrium case with / ~ 6(1), I a ~ T 3 and ~ T 2 , the elastic 
collision rate scales as ~ a 2 s T and the inelastic rate scales also as ~ a 2 s T so they 
are at parametrically the same order; 

3) In the glasmadike overpopulated case with / ~ 6(l/a s ), I a ~ Q^/a 2 . and 
It ~ Q 2 /a s , the elastic collision rate scales as ~ Q s and the inelastic rate scales 
also as ~ Q s so again they are at parametrically the same order. 

We therefore see that the effect of the inelastic collision is parametrically as 
important as the elastic one both near and far from equilibrium, and including 
the inelastic collision qualitatively changes the ultimate equilibrium solution. It 
is clear that even starting from highly overpopulated initial condition, eventu- 
ally there will be no chemical potential nor any condensate in the final thermal 
distribution with the presence of inelastic collision. However, the very impor- 
tant question that has not been understood, is how the inelastic collision will 
affect the transient dynamical off-equilibrium condensation driven by the pure 
elastic evolution starting from initial high overpopulation. Will the system still 
reach the onset of such condensation? Will the inelastic collision speed up, de- 
lay, or completely eliminate such an onset? We will address these questions by 
numerically solving the above kinetic equations. 



3 Numerical Study of the Kinetic Evolution 

In this section, we numerically solve the kinetic equation (18), starting with a 
Glasma-type initial condition as follows: 

f(p,t = 0) = f 6(Q s -p) (22) 

with Q s the saturation scale. We use Q s as unit for all momenta/enegy etc and 
use l/Qs as unit for time. As studied in [16], with a given initial occupation /o, 
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the overpopulation parameter is rte~ 3//4 = f^ 4 2 5 / 4 /3-7r 1 /2 and when f > f§ ~ 
0.154, the system is overpopulated as compared with the Stefan-Boltzmann 
limit and the system will reach onset of condensation when there is only elastic 
collision. For simplicity we fix the initial occupation /o = 1 which is in the 
overpopulated regime. Note that for the the constant t;a 2 a can be absorbed 
by a redefinition of time variable, t — > (£,a 2 )t which we will use from now 
on. We therefore are left with three parameters, R, D, and z c . The inelastic 
contribution will increase with increasing R and a while the elastic will increase 
with increasing D. We will study the effect of inelastic collision by comparison 
with pure elastic case (R = versus R ^ 0) and by varying the strength of 
the inelastic contribution. In passing let us also mention that we have studied 
the case with only inelastic collision, and in that case indeed the system simply 
thermalizes into equilibrium Bose-Einstein distribution with temperature to be 
the supposed value required by energy conservation and without any chemical 
potential. 

3.1 Prom overpopulation toward onset of condensation 

Let us first briefly summarize the kinetic evolution from overpopulation toward 
onset of condensation in the purely elastic case, as reported in [16, 31, 32]. The 
elastic kernel can be re- written in terms of flux as 

C 2 « 2 - ~^d p { P 2 S{p)), S( P ) = - [I a d p f p + I b f p (l + f p )} (23) 

The strong overpopulation leads to a particle flux cascade toward the infrared 
regime. Analysis of the small p regime shows that it will quickly develop a local 
thermal form f*(jp) = \ / [eS p ~ ^ ^ T — 1] with T* — I a /h, and the incoming 
flux will drive the (negative) \i* to eventually vanish and reach the onset of a 
dynamical condensation. This picture is numerically verified in great details in 
[16]. Our discussion of the onset of condensation will stay in this picture (as 
the elastic term is still present and its flux drives the small p behavior) and we 
will study how the inelastic process modifies it. 

Starting with the overpopulated initial condition (22) we have numerically 
evolved the kinetic equations (18), (20), and (21) with given set of parameters. 
Shown in Fig. 2 is the solution with R = l. In Fig. 2 left we show the distribution 
function f p at different time moments, and one can see that even with the 
presence of inelastic term, the small p part of the distribution is quickly filled 
up and becomes a local thermal form f*(p) = l/[e^ p_A1 ^ T — 1] despite that 
the distribution in the wide range of (bigger) momentum region is still far from 
equilibrium shape, and the small p part becomes steeper and steeper with time 
(meaning decreasing |/^*|). In Fig. 2 right we show the corresponding flux S(p) 
from the elastic kernel. Again the flux behaves very similarly to the purely 
elastic case: one see a linear behavior at small p, S cx — p and eventually upon 
onset of condensation the flux diverges (see the blue curve near p = 0) . 

To get an intuitive idea of the contribution of the inelastic kernel, we plot 
the Ci^2 and p 2 C\^2 in Fig. 3. One can see that the kernel is large and 
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p/a p/a 



Figure 2: (Left) The distribution function f{p) at different time moments dur- 
ing evolution; (Right) The flux S(p) defined in elastic kernel at different time 
moments during evolution. 

positive at small p, small and positive at large p, while negative at intermediate 
p. This could be qualitatively understood: significant number of particles with 
intermediate momenta will merge toward high momenta and split toward low 
momenta that will fill up UV and IR region while decrease the occupation at 
intermediate momenta. We also notice that upon onset (the blue dashed curve 
in the right panel) the inelastic kernel near p = shows a divergent behavior in 
consistency with the elastic flux behavior. 

One can directly examine the locally determined T* and /x* (see [16] for 
details) at each time moment during the evolution: these results are shown in 
Fig. 4. Here we also compare the results for different strength of the inelastic 
collision R = 0,0.1,1,10 (noting that the R — case corresponds to purely 
elastic collision). In all cases we can see that the local "chemical potential" /i* 
decreases rather rapidly toward zero. We also show the distribution f(p) at 
the smallest grid point in our calculation p = 0.005Q S as a function of time in 
Fig. 5, which shows very rapid increase of the occupation in consistence with 
the vanishing of (i* . What is most striking is that with increase values of R this 
evolution toward the onset of condensation fi* — > becomes faster and faster. 
The R = 1 case is already much faster than the purely elastic case. This is 
to say, contrary to expectation that the inelastic process may "kill" the strong 
overpopulation quickly, the existence of inelastic collision actually speeds up 
significantly the process of populating the infrared regime and building up a 
local thermal form with vanishing fj,* , which when combined with the structure 
of elastic kernel will then lead to the onset of condensation. 

3.2 Small p analysis of the inelastic kernel 

To understand better the influence of inelastic collision on the small p region, 
let us examine the kernel (21) for p — > limit before the onset of condensation. 
Provided that /o = f(p — 0) < oo and that its derivatives with respect to p at 
p = is also finite, we can have the expansion f(p — > 0) « fo + f^p + .... If 
we look at a small enough p regime, then all the involved momenta (p, zp, (1 — 
z)p,p/ (1 — z),pz/ (1 — z)) in the kernel (21) can be considered small and we can 
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Figure 3: The inelastic kernel Ci^ip) (Left) and p 2 Ci^2(p) (Right) at different 
time moments during evolution. 




Figure 4: The local thermal form parameters T* (Left) and /i* (Right) as 
functions of time for different values of parameter R. 
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Figure 5: The occupation at the smallest grid point f(p = 0.005Q S ) as a function 
of time for different values of parameter R. 

use the expansion for them. This leads to 

Cft a (p->0)->fl£ U/o(l + /o)+Ai^(l + 2/o)p + 6(p 2 )l , (24) 

where we have introduced the constants 

1 1 z c {\\z 2 c - 27z c + 18) 



A = In 



1 - % 6 (1 - z c f 



1 1 z e (25z c 3 - 88z c 2 + 108z c - 48) 

Al ln —-12 (HT^ ' (25) 

All these A's are positive for < z c < 1. Clearly for sufficiently small p the 
leading term in the inelastic kernel ~ /o(l+/o)A) is always positive and becomes 
bigger and bigger with increasing /o (which is a kind of "self-amplification"). 
This will tend to increase the particle number near p = very rapidly and the 
effect becomes stronger with increasing values of R, which explains the behavior 
seen in Fig. 5. 

Physically this behavior may be understood in two ways. First note that the 
inelastic kernel has its fixed point to be \ j(ePl T — 1) which at small p is ~ l/p 
so as long as f(p = 0) is finite yet the inelastic kernel will try to fill it up toward 
l/p. Second, this is also related to the Boson nature: if all involved particles are 
from small p, then the merging rate is like ~ f§ (1 + /o) while the splitting rate 
is like ~ /o(l + /o) 2 so the splitting "wins" due to Bose enhancement for the 
final state particles and it increases particle number at small p. To conclude, 
the inelastic kernel contribution is always positive at very small p and it will 
catalyze and speed up the onset of a Bose condensation (which itself is driven 
by the elastic term at fj,* 0). 



3.3 Change of particle number from inelastic kernel 

While the inelastic kernel always increases the occupation at sufficiently small 
p, it may still decrease the the total particle number. Indeed as shown in Fig. 6 
(left panel), the total particle number decreases when R > 0, and it decreases 
more rapidly for larger R. 
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Figure 6: (Left) The particle number density as a function of time for different 
values of parameter R. (Right) The entropy density as a function of time for 
R = 1 and R = 0. 

To understand the change of particle number n = J p f(p), one can integrate 
the two sides of the kinetic equation (18) and obtain 

d t n(t) = R^j b J dp? 2 J Q YZ~ z [f P 9(i- z)p9zp 9pf{l-z)pfzp\ 

= R 2^h I ^ L ^ P + ^ (1 ~ 2)p + ^ zp ~ fQ--*)pf*p[ -( 26 ) 

From the above one can see the for the region z — > the leading order in the 
z-integrand becomes ~ /p(l + /p) and the contribution is positive, i.e. increasing 
particle number. For general z, the z-integrand can be rewritten as 

dMt) = R 2^r b j dpp2 r T^~ z [/p(i + fp) ~ {f(i - z)p ~ fp){fzp ~ • (27) 

We see that for not too small z, the momenta zp, (1 — z)p become well separated 
from p and the second term in the above integrand becomes important and its 
contribution is negative which decreases the particle number. 

In Fig. 6 (right panel) we also show the entropy density as a function of time 
and compare the case with R = 1 and the purely elastic case with R = 0. One 
can see that with inelastic collision included, the entropy density increases much 
faster. That is, the inelastic process tends to accelerate the thermalization. 



3.4 Dependence on the parameter z c 

Finally we study the dependence on the parameter z c which is the kinematic cut 
to make sure the validity of the approximations used for the matrix element. 
Generally speaking, with larger z c we include more effects from the inelastic 
process. To see how the results depend on z c , we fix other parameters and 
compare the results for different choices of z c . In Fig. 7, we show the local 
thermal form parameters T* (left panel) and fx* (right panel) as functions of 
time for z c = 0.2, 0.5, 0.8 respectively. In Fig. 8, we show the occupation at 
the smallest grid point (left panel) and the total particle number (right panel) 
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Figure 7: The local thermal form parameters T* (Left) and //* (Right) as func- 
tions of time for z c = 0.2, 0.5, 0.8 respectively. 
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Figure 8: The occupation at the smallest grid point f(p = 0.005Q S ) (Left) and 
the particle number density (Right) as functions of time for z c = 0.2, 0.5, 0.8 
respectively 



14 



as functions of time for z c = 0.2,0.5,0.8 respectively. From the plots we can 
see that indeed with larger z c the fo increases faster and /i* vanishes faster as 
expected for stronger inelastic effect. For the particle number, the case with 
z c = 0.2 actually has n increasing with time, which can be understood fro m 
the analysis in the previous subsection. The particle number in both z c = 0.5 
and z c = 0.8 cases drops with time and does so faster for larger z c . In passing 
let us mention that we have also studied the dependence on the parameter D: 
basically increasing D will enhance the effect from elastic collision and also 
speed up the thermalization in general, as well as reach onset of condensation 
at earlier times compared with D = case. 

4 Conclusion 

In summary, we have studied the kinetic evolution of a highly overpopulated 
system starting from Glasma-type initial condition with the presence of both 
elastic and inelastic collisions. Using the Gunion-Bertsch formula for the 2f>3 
matrix element, we have derived the inelastic collision kernel under the collincar 
and small angle approximations. Putting together the inelastic kernel together 
with the previously obtained elastic kernel, we have then numerically solved 
the kinetic evolution for varied choices of parameters. Our main finding is that 
the inelastic process has two effects: globally changing (mostly reducing) the 
total particle number, while locally at small p always filling up the infrared 
regime extremely quickly. The latter effect is shown both from numerics and 
by analytic analysis. This effect significantly speeds up the emergence of local 
thermal form near p = and the vanishing of local "chemical potential" fx* 
as previously found in the purely elastic collision case to lead to the onset of 
dynamical Bose condensation. Therefore in our present approach of including 
the inelastic scattering, we conclude that, contrary to some previously discussed 
expectations about the role of number non-conserving processes, the inelastic 
collision actually helps to build up the local "critical form" ~ 1/p much faster 
and catalyzes the onset of condensation in the overpopulated Glasma. 
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Figure 9: A typical Fcynman diagram for 

Appendix 

A gg — > ggg matrix element square 

The invariant gg —> ggg matrix element square summed over all final states 
and also summed over all initial states is computed by considering 25 different 
diagrams [40]. We quote it here: 

|Af 2 ^ 3 | 2 = g 6 N*N g ^ [(12345) + (12354) + (12435) + (12453) + (12534) 

+(12543) + (13245) + (13254) + (13524) + (14235) + (14325)], 

(28) 

where N c = 3 and N g = 2(N 2 - 1), 

N = (12) 4 + (13) 4 + (14) 4 + (15) 4 + (23) 4 + (24) 4 

+ (25) 4 + (34) 4 + (35) 4 + (45) 4 , 
V = (12)(13)(14)(15)(23)(24)(25)(34)(35)(45), 
(ijklm) = (ij)(jk)(kl)(lm)(mi), 

(ij) EE ki-kj. (29) 

|M 2 ^.3| 2 is completely symmetry in fej, i = 1 — 5, so let's take k\ and fc 2 as the 
hard momenta in the entrance channel, k% and k 4 as the hard momenta in the 
exit channel, and k§ as the emitted soft gluon. We also denote the exchanging 
momentum as q = fc 2 — k 4 . A typical Feynman diagram illuminating this set-up 
is shown in Fig. 9. 

Define the Mandelstam variables as [41, 42]: 

a = (fci + k 2 f = 2(12), t=Qs 1 - k 3 ) 2 = -2(13), u = (h - h) 2 = -2(14), 
s' = {k 3 + k 4 ) 2 = 2(34), t' = {k 2 - k 4 ) 2 = -2(24), v! = (k 2 - k 3 ) 2 = -2(23). 

In addition, the following relations involving k$ hold: 

= B + f + « = s + t' + u'^ = s + t' + u^ = s + t + u'_ 
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In terms of the Mandelstam variables, |M2_>3| 2 can be written as 



\M; 



2^3 I 



32g 6 N*N g Af 



1 



1 



1 



s'(s + u + t)(s + u' + V) \W uu' 
1 /ll 



where 
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+t>)\ 


i'(s4 
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f« 


f*)C 


t(s + 


u' 


+*)(*- 
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+- u' 


+ t')( 


u'(s - 


V u 


+ t)(s 
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+ u' 


+ t)\ 


u(s 4 
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+■ 


+ u> 


+ t')( 


[s 4 + 


t A 


-f u 4 + 




ht' 4 + 


s + u 


' + 


t') 4 + 


{s + 


t' + u) 



(30) 



.Ml 



Assume q = k 2 — k± is small as well as k§. Then t, t' are small, while s, s', u, v! 



are large, and — u 



s' —> s. We do an expansion for \M\ and keep 



only leading and sub-leading order terms in q and k$. First, 

s' = (ki + k 2 - fc 5 ) 2 = s - 2(fci + fc 2 ) • k 5 + 0{k 2 ), 

u = -s + 2k\ ■ k 5 + 0{q 2 ), 

u' = -s + 2k 2 -k 5 + 0{q-k 5 ) 2 , 

t = (q-h) 2 , 



Second, 



fe ■ (h + q) = => k 2 ■ q = -k 2 ■ k 4 = -(k 2 - k 4 ) 2 = -q 2 , 
t-t' = -2q ■ k 5 => -2q ■ k 2 - 2ki ■ (q - k 5 ) = 2k 2 ■ fc 4 - 2ki ■ k 3 = t - t' = -2q ■ k 5 



h ■ (k 5 - q) 



-q-h 
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Then Af and |M2_).3| 2 can be simplified as 



Af = _- s J (fc 2 .fc 5 + fc 1 -fc 5 )+0(fc 5 2 ,g i ,fc 5 -g), (31) 



1 1 



s'(s + u + t)(s + u' + t') s(s + u' + t)(s +u + f) 



1 1 



u'{s + u + t)(s + u' + t) u(s + u + t')(s + u' + V) 
qo 6 M 3 M M s + (k 1 +k 2 )-k 5 2 
32gN ^ ^k 1 .k 5 )(k 2 .k 5 ) + °^ q) 

^ 6 ^ g ^ s - 3( M fc2) ,;t 5 +Q(fc^) s 



o 




'tt' (h ■ k 5 )(k 2 ■ k 5 ) 

= ^?N g <** ' % 2fcl ; - + ^ ' fc5 + 0{k 5l qf. (32) 
q 2 (q-k 5 ) 2 (h ■ k 5 ){k 2 ■ k 5 ) 

In the center of mass frame of ki and fc 2 , it goes to 

|M 2 »,f * |Mqb| 2 fl - 3^ fl + *^*) fl + te ',:f;^ ) ; k5ll)2 

|M GB | 2 (^l-3 (33) 

where q|| = (q • vi)vi and 

\M GB \ 2 = 32 g 6 7V e 3 iV "I (34) 

is the Gunion-Bcrtsch formula. We used the fact that, in any frame, q$ = 
(v 2 ■ q) 2 + 0{q\q a ) and (q - fc 50 ) 2 = (vi • (q - k 5 )) 2 + 0(q± - k 51 _) 2 (q - fc 50 ) 
for soft q and fcg. Thus, Gunion-Bertsch formula is the leading order amplitude 
in soft q and k expansion; and Eq. (32) is the amplitude including both the 
leading (Gunion-Bcrtsch) and the sublcading order terms. Higher order terms 
can also be obtained, but we won't use them. Note that one can naively boosts 
the Gunion-Bertsch formula from the center of mass frame to a general frame 
by using the replacements k 2 ± — > 4(ki ■ k 5 )(k 2 ■ k$)/s, q x — > 4(fci • ki){k 2 ■ k^)/s, 
and (q_i_ — k 5X ) 2 — > A{k\ ■ k 3 )(k 2 ■ k 3 )/s. The result coincides with Eq. (32) at 
leading order but not at next to leading order in soft q and k$ expansion: 

\M GB \ 2 ^l6g 6 N^N g - r~ul i Ui ^TT—i Ui FTn ( 35 ) 

(ki ■ h)(k 2 ■ k 3 )(ki ■ k 4 ){k 2 ■ fc 4 )(fci • h)(k 2 • fc 5 ) 

B The collision kernel £2^3 at collinear approx- 
imation 

The collision kernel C 2 ->3 has a very complicated structure, in this section, we 
simplify it by taking the collinear approximation, i.e., Vfc ~ vi or ~ v p . 
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We thus rewrite the collision kernels C 2 ^ 3 and C 2 ^ 3 as (We denote 1 + /; 
by 9i) 

x{2nfS\l - k - p 3 )\M^ 23k \ 2 (g pgi f 2 f 3 f k - f p fig 2 939k), 

**■ " £ L, H / (0 / t< 2 *>'^ + « - '^'^ - « - » > 

x (2tt) 4 5 4 (Z + fc - P3)|M 2 h 3^ lpfc | 2 (.g p5l5fe / 2 / 3 - fphhgnga), (36) 

where we introduced two auxiliary integrations over / and g, the kinematics is 
shown in Fig. 1. 

First, it is easy to show that, under small angle approximation, 

,/2 



q Q * q-V!-^ 



(qx - kj ^ 
,//2 



^ + (q-k).v p + V ^ oF ^ , (37) 



2E p 



h ^ £i + Sfc - k • v ; 



£ P + 9o, (38) 



. q_L qi . (q-v p ) 2 q-Vp 

v < " Vp + k m Vp + ^r v - (39) 



^p ^J^p J-^p J^p 

where qj_ = q q • v p Vp, kj_ = k k • v p v p , q^_ = q q • vivi, and k"j_ = 

k k ■ V;V;. 

Second, if k < q, then, under collinear approximation, is either nearly 
parallel to v p or vi. For v fc ~ vi, q ■ k = |k|(q - q ■ v fe ) ~ \k\[q - q ■ vi + q ■ 
(vi-V fc )] ~ -|k|q'l/^ + |k||q|0 lfc < q 2 ; For v fe ~ v p , q-k = |k|(|k| -k- v,) ~ 
|k| 2 ^. p + \k\cfi_/(2E p ) + |k||q|0 fep < q 2 . Thus in the small angle plus collinear 
approximation, if k < q, we can approximate the amplitude \M1 p ^ 2 ^ k \ 2 as 



For |-^23->ipfc| 2 ' U P t° 9 2 order, we can neglect q ■ k in the denominator and g -p 
and q • pi in the numerator, 



\M b 23 ^ lkp \ 2 = Mg e N?N 9 2I [ t^ q) ; {P X k uif ^ n 
1 23 ^ Upl y c 9 g 2 (g + fc) 2 [(p + fc- ( 7)-A : ][(pi+ 9 )-fc] 



9 (g 2 ) 2 [(p-9)-*][(Pi + g)-A] 



9 (q 2 ) 2 (p-k)( Pl -k) 



\M^ 23k \ 2 . (41) 
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Third, if k > q, under the collinear approximation, for Vfe ~ Vi, q ■ k ~ 
-|k|q'l/Ei + |k||q|6»i fe < q 2 because 6" lfe < Iq'U/^i! Fot v fe ~ v p , g • fc - 
|k| 2 (9 2 p + |k|q2_/(2£;p) + |k||q|0 fe p < g 2 because 6 kp < |pj/£ p . Thus the 
collinear approximation can simplify the amplitude even when k > q: 



|M&_> a3fc | a - 64. 9 6 iV r 3 7V 9 



Fourth, the whole kinematic (phase) space can be separated into two parts: 
one for k < q and the other for q < k. We would expect that at kinematic region 
with k < q, the 2 — > 3 process may be regarded as a 2 —> 2 "hard" process with 
one additional "soft" gluon emitted or absorbed by one of the "hard" gluons; 
with q < k the 2—^3 process can be regarded as an effective 1 — >• 2 process 
with a spectator gluon joined to make the effective 1 — > 2 amplitude nonzero 
(the amplitude of f —> 2 process is zero for massless gluons). Thus we separate 
the collision kernel as 

C2^3 = C^ 3 + C 2 ^ 3 = C^ 3 + C^ 3 , (43) 

with 

x |M 1 ^ 23fc | 2 [(2^) 4 ( 5 4 (/ - k - P z){g pgi hhfk - f P .fi92939k) 
+(2tt) 4 (5 4 (/ + k - Pa)(g P gigkhfa - Uhh92gz)\, 

x [{2vf5\l -k- Pa)\M^ k \\g vgi hhfk - f P fi929z9k) 
+ (2tt) 4 S 4 (1 + k-p 3 )\M^ lkp \ 2 (g p g igk f 2 f 3 - f P hhg 2 ga)]. 
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B.l Simplifying C 



> 

2^3 



Expand the integrand of in terms of k and keep the leading order terms: 
» 1 f di 1 f f q>k 1 f dl °,o ^8,4, , ns i, 



«- "55 y p# l l ^ y ^ <2 " 5 (p + « - m <" - 9 - b) 

x|A/i%23fe| 2 [(27rWo - l k l - I 1 - Ul)(9pgihfif* ~ Uh929l9k) 
+ (2^)<5(Z + |k| - |1 + kDGfrfftffc/a/, - f P hh929i)\ 

' dh {2,f5% + q -l)8\ Pl - q - P2 ) 



N g J (27r) 4 A 2 ;A 2S p y 2^ 
x|A^23 fe | 2 [(2T)^(^o - |l|)(fl P fli/ 2 /j/fc - f P fi929Wk) 
+ (2n)6(l Q - \l\)(g p gi9kf2fi - f P hfk929l)] 
d 4 n f r q>k 1 

-(2tt) 4 ( 5 4 (p + 9 - 0(2tt) 4 5 4 (pi - ? - pa) 



AW (2n)*JwJk 2E . 



. 2 2g 2 N c {p-p x ) i 



x\M lp ^ 2l \' " , fc) T ^ (1 + Vk){9p9lhfl - f P fx929i), (44) 



where in the last equality Z M = (|1|, 1) is on-shell, and the 2 — >• 2 amplitude is 

\M lp _> 2l \ 2 = 32 ff 4 iV c 2 7V 3 %^. (45) 



We have written the 2 — > 3 amplitude in a form of a 2 — > 2 amplitude times a 
1 — > 2 splitting function. Indeed, if for example k is nearly collinear to p, 



2g 2 N c ( P - Pl ) 2g 2 N c \v\ 



(p-fcXpx-fc) (p-fc)|k| 
2g 2 
(p + *) 



2^(^) (46) 



with Pgg(z) = 2Ca/z being the standard unregularized g — » splitting func- 
tion at z -> limit where z = E k /E p [43, 44, 45]. 

In the collinear approximation, if Vk — vi, then p ■ pi/k ■ p ps |pi|/|k|; or if 
Vk ~ v p , then p ■ pi/k ■ pi ps |p|/|k|. We have 



ly g J\2l LEj V Jk< P1 ~p 2 |K 



2g 2 N c 



1 + 2/ fc 1 + 2/ fc 



1 — Vfc • Vi 1 - Vfc ■ Vp 

x(27r) 4 ^ 4 (p + pi -p 2 -l)(g P gif 2 fi - fpfmgi)- (47) 

This is essentially a.p + pi p 2 + l collision kernel with an inner 1 — > 2 splitting 
function. Let 



2%) 



2g 2 N c 
k< q |k| 2 



1 + 2f k 1 + 2/ fc 



1 - Vfc • vi 1 - v fc • v p 



(48) 
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For isotropic distribution, 

d 3 k 2g 2 N c (l + 2f k ) 



V{q) = 2 / 

Jk 



k<g (2») 3 2E* M 2 (l - cosfl) 



(2tt) 2 J 1 1 |k| J 1- cos 6 
a 2 N f' q ' I + 9 f, f 1 

>9 iVc / dlkli-^F / (49) 



(2tt)2 y 1 1 | k | j_, i 



Thus, 



<£+3 » / ^r|M lp _ >2; | 2 P(pi-p 2 )(2^) 4 5 4 (p + p 1 -p 2 -/) 

x(g P gif2fi - fphgigi)- (so) 

There are two types of infrared divergence in T>(q). (1) The logarithmic diver- 
gence J dx/(l — x) ~ J e?#/6> ~ ln(l/0 m ) with # r „ the minimal angle between 

k and p. m arises completely due to interaction, so m ~ g. Thus, J_ 1 da;/ (1 — 
a;) ~ ln(l/<7) in both glasma and nearly thermal equilibrium. (2) Near thermal 
equilibrium, f k ~ T/w*. thus /*(dfc/fc)(l + 2/ fe ) ~ 2Tj^(dkk/(k 2 + m^) 3 / 2 ~ 
T(l/m oc — 1/mjj) where we use moo to denote the mass of the emitted or ab- 
sorbed ultrasoft gluon fc and to denote the mass of the exchanged gluon 
q. Near equilibrium, both m,oo and mo are of order gT but can have different 
prefactors, we find T>(q) ~ gln(l/g). In the initial glasma, mo ~ moo ~ Qs 
and h ~ l/a s , thus f*(dk/k)(l + 2f k ) ~ (2/a.) In(g/moo) - 1/a,. Thus 
V{q) ~ ln(l/g). As the glasma evolves, ~ As/(a a k), if m_o ~ moo ~ VAAg, 
thus / g (d fc/fc)(l + 2/fe) ~ (l/a s )(A s /moo) ~ {l/a a )yfK^JK. Thus, we find 
P(g) ~ v^7Aln(l/ 5 ). 

In any case, we can conclude that the ratio of ultrasoft gluon emission and 
absorbtion 2 — !► 3 processes over the elastic 2 — > 2 processes is either ln(l/g) 
order or g]n(l/g) order. 



B.2 Simplifying C 



< 

2^3 



Expand the distribution functions in C^i->3 m terms of q and keep the leading 
order terms: 

x [(2^) 4 <5 4 (p + q - k - P3)\M^ 23k \ 2 (g p f 3 f k - f p g 3 g k ) 

+(2tt) 4 5 4 (p + q + k - P3)\M^ lkp \ 2 (g P gkf 3 - f P .fkg 3 )} 

IIW^L Ik" U&ZWto - Q ■ vOIM^fe l^r 
x [(2^) 4 <5 4 (p + q-k- P 3 )\M^ 23k \ 2 (g p f 3 f k - f p g 3 g k ) 
+(2tt) 4 5 4 (p + g + fc - P 3 )\M^ lkp \ 2 (g p g k f 3 - f P f k g 3 )}, 



N, 
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where hi = fig\. This is basically a 1 -> 2 collision kernel with the 2 — 2 
processes playing a role of opening a finite phase space for 1 — > 2 process. 
Because q is small, p, k and p$ are nearly collinear (there is no phase space at 
the collinear region k \\ p\). In the collinear approximation, we have 



p a + q Q - k - p 30 
Po + qo + h) - P30 



q - |k| - (q - k) • v p « q - q • v p , 
<7o + |k| - (q + k) ■ v p w q Q - q • v p . 



The amplitudes are 



lp— J-23A; I 



l-^23-!-lfcp| 2 



64g 6 N*N, 



(P-Pi) 2 



s (g2)2|k|2(l_ Vp .v fc )' 
(P-Pi? 



Thus, 



L '2->3 



AT,, 



d 4 g 



64g 6 N^N g ^- 



q<k 



k| (q2)2| k |2 (1 _ Vp . Vfc) - 



(51) 

(52) 



1 Jk 



(2E P )*2E : 



(2Tr)8(q - q • v 1 )(2ir)5(q - q • v p )/ii 



x [iM^agfcl^pp/p-fe/fc - f p g p -k9k) + \M% 3 ^ lkp \ 2 (g p gkf P +k - f P hg P +k)\ ■ 

(53) 

In the following we denote g = |q|,p = |p|,Pi = |pi|, k = |k|. Let qo = x<7 and let 
v p = (1, 0, 0), vi = (cos 0i, sin 0\, 0), and v g = (sin# g cos , sin # g sin 9 , cos q ) . 

5(x — v q ■ v\)6(x — v g ■ Vp) = S[x — sin6* g cos(#i — <p q )]5(x — svn0 q cos0 9 ) 
1 S(x — sin # g cos 4> q ) 



s\\\0 q | sin0 g — sin(<^ g — 



1 



1 



sin6» 9 2sin(6»i/2) 



S 4>, 



Six — sin n cos ■ 



+8 \4> q — ^ — 7r J 5 ( x + sin q cos -J- 



Thus, 



| ^fi^Tr)^ - q ■ Vi)(2n)S(q - q • v p )|A^ 23fc 



1287r ff 6 iV c 3 iVg 



q 5 k 2 
(l-vp-vi) 2 



c/.r 



d0 g sin g 



<5(.t - v, • vi)S(x - v 9 • v F 



(l-x 2 ) 2 (l-vp.v fc ) 

(ppi) 2 (l-vp-vx) 2 



= l28irg 6 N^N g - 



= l28Tr 2 g 6 N^N g - 



q 5 k 2 1 - v p • Vfc 

(ppi) 2 3 - Vp ■ vi 
<j 5 /c 2 1 - Vp ■ va_. 



d0r 



sin^(l - sin 2 q cos 2 %) 2 



(54) 
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Furthermore, for isotropic distributions, we have: 



128n 2 g 6 N%N g 



dqq 2 



q<k 



1 Jk 



1 (ppi) 2 3 - v p ■ vi 
(2E p ) 2 2E 1 q*k 2 l-v p -v fc 



{gpfp-kfk - f P g P -k9k) h 3 — {g P 9kfp+k - f P fkg P +k) 



3g 6 N 3 
16tt 5 



dpxp\h\ 



dx 



dq 



x Jo F 



dk 

T 



(p _|_ fc)3 

(g P f P -kfk - f P g P ~kgk) H 3 — (g P gkf P +k - f P fkg P +k) 



,(55) 



where the upper limit of the integration over fc for the first two terms should be 
cut at p. When k is small the integrand over k goes like (1 + 2f p )f'. It is finite, 
so we can put the lower limit of J dk as 0. In the first two terms, let k = zp 
with z being the momentum fraction of the emitted gluon; in the last two terms 
let k = z(p + k) with z being the momentum fraction of the absorbed gluon. 
Then we have 



t -'2-J>3 



1 



dz 
z 



{i-zY 



[g P f(i-z) P fz P — f P g(i- z ) P gz P ] 

[9p9zp/{i-z)f P /(i-z) - fpfz P /(i-z)9p/(i-z))] \, (56) 



where we introduce the momentum fraction cut z c < 1 to characterizing that 
k being a small fraction of the total momentum in this effective 1 — > 2 process 
and the prcfactor R is given by 



R 



127V 3 1 



dx 



t J-il-x 



dq 

1 



0(1). 



(57) 



There are two kinds of infrared divergences in C^^' (1) The logarithmic diver- 
gence: J dx/(l — x) ~ ln(l/#). (2) The quadratic divergence J dq/q 3 ~ l/m 2 D . 
Noticing that C2->2 is of order a 2 \i\{l/g) near equilibrium and (Q s /a s ) ln(l/<7) 
in initial glasma state, we find that C^_ >3 /C2^2 is of order O(l) in both equilib- 
rium and initial glasma states. 

Now we show that the collision kernel C^_j. 3 conserves energy, i.e, 

So <h>iK 2 

dpp 3 C 2 < ^ 3 [f P } 



0: 



dpp 3 
1 



~/ [g P f(i-z) P fz P — f P 9(i-z)p9zp] 
4 [g P gz P /(i-z)f P /(i-z) - f P fzp/(i-z)9p/(i-z))] 



In the last two terms, change the variable p — > (1 — z)p, one finds that the first 
two terms cancel the last two terms so that L dpp 3 C^ 3 [/ p ] = 0. 
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